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ABSTRACT 

We present the first scattered light images of debris disks around a K star (HD 53143) and an F star 
(HD 139664) using the coronagraphic mode of the Advanced Camera for Surveys (ACS) aboard the 
Hubble Space Telescope (HST). With ages 0.3 - 1 Gyr, these are among the oldest optically detected 
debris disks. HD 53143, viewed ~45° from edge-on, does not show radial variation in disk structure 
and has width >55 AU. HD 139664 is seen close to edge-on and has belt-like morphology with a dust 
peak 83 AU from the star and a distinct outer boundary at 109 AU. We discuss evidence for significant 
diversity in the radial architecture of debris disks that appears unconnected to stellar spectral type or 
age. HD 139664 and possibly the solar system belong in a category of narrow belts 20—30 AU wide. 
HD 53143 represents a class of wide disk architecture with characteristic width >50 AU. 

Subject headings: stars: individual(HD 53143, HD 139664) - circumstellar matter 



1. INTRODUCTION 

The configuration of our solar system is perhaps the 
most significant starting point for our understanding of 
planet formation. Therefore a fundamental question is 
whether or not the architecture of our solar system is 
common relative to other planetary systems. One point 
of comparison is the structure of our Kuiper Belt relative 
to other systems, which are typically seen as debris disks 
in scattered light or thermal emission. In scattered light, 
some debris disks, such as (3 Pic and AU Mic, have central 
holes, but are radially extended to hundr eds of AU radii 
l|Smith fc Terrilelll984t Ikalas et alJ l2b041. Other debris 
disks, such as HR 4796A and Fomalhaut consist of rel- 
ativel y narrow rings with sharp inner and outer bound- 
aries (IHSEnildiZilin] EilaZiLaniinnii)- How- 
ever, a narrow-belt architecture has not previously been 
detected in scattered light among stars similar in spectral 
type and age to the Sun. 

HD 53143 (K1V) and HD 139664 (F5V) are two 
stars ~18 parsec from the Sun known to have cir- 
cumstellar dust due to excess thermal emission at far- 
infrared wavelengths (Aumann 1985; Stencel & Back- 
man 1991; Table 1). Various in dicators place the 
age of HP 53143 at 1.0 ± 0.2 Gyr iPecine et alJl2000l 
iSong et all 1200(1 iNordstrom et alJ l2004|h whereas HP 
139664 may be a younger system with age 0.3^2 Gyr 
(iLachaume et al.lll999t [Montes et~aTll200H iMallik : et all 
120031 INordstrom et al]l2004l) . Forthese two stars the 
infrared excess corresponds to a dust mass 3-10 times 
smaller than t hat of ~10 Myr-o ld systems such as AU 
Mic and (3 Pic (|Kalas et al.ll2004l) . Pirect imaging of de- 
bris disks with masses this small is observationally chal- 
lenging, but it is now feasible using the optical corona- 
graph in ACS. 

2. OBSERVATIONS & DATA ANALYSIS 

We utilized the HST ACS High Resolution Camera 
(HRC) with a 1.8" diameter occulting spot to artificially 
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eclipse each star (Table 1; Fig. 1). Five F stars were 
observed in consecutive orbits, as were five K stars, in 
order to minimize differences in the point spread func- 
tion (PSF) due to telescope thermal variations. Each 
PSF was subtracted using the four other stars in each 
set of observations. The relative intensity scaling and 
registration between images was iteratively adjusted un- 
til the residual image showed a mean radial profile equal 
to zero intensity. 

After the excess nebulosity was detected around HP 
53143 and HP 139664, we determined that no surface 
brightness asymmetries were detected between each side 
of each disk. To improve the signal-to-noise, we mir- 
ror averaged the data (Figs. 2 & 3). Mirror averaging 
splits the image into two halves along the axis that is 
perpendicular to the disk midplane and bisects the disk. 
One side is transposed onto the other side and the data 
are then averaged. Asymmetries due to scattering phas e 
function effects will be coadded l)Kalas fc Jewittl il996). 
In effect mirror averaging doubles the integration time 
spent on the circumstellar disk given that the broad fea- 
tures between each side are symmetric. As a test, we 
also subtracted the two disk halves from each other and 
confirmed that the assumption of symmetry is valid. 

3. RESULTS 

The two disks have different morphologies due to differ- 
ent inclinations and intrinsically different architectures. 
To quantify the viewing geometries and structural prop- 
erties of the disks, we produce a series of simulated 
scattered-light disks that explore the parameters of in- 
clination to the line of sight, inner and outer disk ra- 
dius, and the ra dial and vertical variati on of dust number 
density (Fig. 2: lKalas fc ,Tewittlll996t) . We reinsert the 
simulated disks across each star in a direction orthog- 
onal to the observed midplanes and select those models 
that most closely resemble the properties of the observed 
disks. Table 1 summarizes our findings. 

The shape of the midplane surface brightness distribu- 
tion differs significantly for each system (Fig. 3). The 
HP 53143 midplane surface brightness decreases mono- 
tonically with projected radius, approximately as r~ 3 , 
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where r is the projected radius. In the simulated disk, 
the radial number density distribution decreases as q^ 1 , 
where q is radius in the disk cylindrical coordinate sys- 
tem. Our solar system's Zodiacal dust complex has a 
comparable dependence of grain number density as a 
function of radius (g~ 134 ; Kelsall et al. 1998), con- 
trolled mainly by the force of Poynting-Robertson (PR) 
drag that causes small grains to spiral into the Sun 
l|Burns et aljfl979^ . Therefore the HD 53143 disk may 
simply represent a population of unseen parent bodies 
that collisionally replenish dust that is redistributed ra- 
dially by PR drag. Our simulations show that the outer 
radius of the observed disk is a sensitivity-limited value 
at approximately 6" radius (110 AU). The inner radius is 
also sensitivity- limited to 3" radius (55 AU). Therefore, 
the debris disk around HD 53143 is at least 55 AU wide. 

Material surrounding HD 139664, on the other hand, 
is confined to a narrow belt, as indicated by a turnover 
in the midplane surface brightness profile between 4.5" 
and 5.5" (79 - 96 AU; Fig. 3). Our disk simula- 
tions show that the peak in the dust distribution oc- 
curs at 83 AU, decreasing as q~ 2 5 from 83 - 109 AU, 
and with a sharp outer truncation at 109 AU (Figs. 
2 & 3). We tested model disks that have outer radii 
>109 AU and found that these disks would have been 
detectable in our data as far as 10" radius (175 AU; 
Fig. 3). Given an absence of significant gas, the belt- 
like nature of HD 139664 is most likely a structure re- 
lated to planet formation. iKenvon fc Bromlevl (2004) 
find that a dust belt with peak surface brightness at 
~80 AU radius forms within a planetesimal disk at age 
400 Myr. The appearance of this belt signals the re- 
cent formation of a >1000 km plane t that gravitat i onally 
stir s planetesimals in its viscinitv. [Liou fc ZooU i|1999f) 
an dM oro-Martin fc Malhotral l|2002fl . on the other hand, 
simulate the concentrations of dust in trans-Neptunian 
space that arise due to trapping in mean motion reso- 
nances. A natural explanation for the belt-like morphol- 
ogy of HD 139664 is that the ~83 AU peak in the dust 
distribution corresponds to either an interior or exterior 
mean motion resonance created by a companion to HD 
139664. Large grains may dominate the belt's peak at 
83 AU, with smaller grains passing quickly through the 
resonance regions due to radiation forces. 

4. DISCUSSION 

Taking a census of eight debris disks resolved in scat- 
tered light and the predicted dust distribution in our 
Kuiper Belt, we observe two basic architectures that are 
not correlated with stellar mass and luminosity, but must 
depend on other environmental factors (Table 2). Debris 
systems are either narrow belts or wide disks. Both types 
have central dust depletions, but the distinguishing char- 
acteristic is the presence or absence of a distinct outer 
edge. HR 4796A, Fomalhaut, HD 139664 and the Sun 
are examples of narrow-belt systems. The belt systems 
appear to have radial widths ranging between 20 and 30 
AU, and the inner edges may begin as close as 25 AU 
(Sun), or as far as 133 AU (Fomalhaut). HD 32297, f3 
Pic, HD 107146, HD 53143 and AU Mic are examples of 
disks with sensitivity-limited outer edges that imply disk 
widths >50 AU. 

The F, G, and K stars are interesting because a priori 
we might expect to find planetary systems similar to our 



own, yet we discover significant diversity in the outer re- 
gions that correspond to Neptune and our Kuiper Belt. 
With age ~1 Gyr, HD 53143 is among the oldest known 
extrasolar debris disks, yet its wide-disk architecture re- 
sembles that of the ~10 Myr old systems of (3 Pic and 
AU Mic. The lingering dust mass (Table 1) throughout 
the system could signal the absence of giant planets that 
otherwise sweep clear the parent bodies (comets and as- 
teroids) responsible for the producing the dust disk. Yet 
the presence of a dust disk out to at least 110 AU ra- 
dius shows that the primordial circumstellar disk prob- 
ably contained the prerequisite mass of gas and dust to 
form giant planets. By contrast, the younger system HD 
139664 has already developed a narrow-belt architecture. 

Narrow-belt architectures for the underlying popula- 
tion of planetesimals may originate from early stochas- 
tic dynamical events, such as a close stellar fly by, that 
strip disk mass and dynamically heat the surviving disk 
<|Ida et al.M200T)t lAdams fc Laughlinll2T)0l . Theoretical 
simulations show that a reduction in disk mass, combined 
with dynamical heating, produces a less stable planetary 
system that is more likely to eject giant planets from 
their formation site to much larger radii, as has been pro- 
posed for the origin of Neptune ijThommes et aUll.9991 
iTsiganis et al.ll2005j) . Therefore, planetesimal belts not 
only evolve into a narrow structures because of exter- 
nal stochastic events, but they may be found at large 
distances from the central star due to the subsequent 
outward migration of interior planets. 

The collisionally replenished dust population will 
spread away from any narrow belt of planetesimals. If the 
scattered light appearance continues to manifest a nar- 
row structure, then both the inner and outer edges are 
probably maintained by other gravitational perturbers 
such as stellar or sub-stellar companions. However, only 
HR 4796 has a known stellar companion that may trun- 
cate the outer radius of the dust belt l)Augereau et alJ 
1999). If there is no confinement mechanism for the outer 
radius, then the architecture will manifest as a wide-disk. 
For example, (3 Pic and AU Mic may have narrow belts of 
plane tesimals ijAugereau et al.l|2001l iStrubbe fc Ghianel 
|2005|) . but the observed dust disk widths extend to hun- 
dreds of AU. Though the predicted structure of our 
Kuiper Belt places the solar system among the narrow 
disk architectures, it is conceiv able that the dust com - 
ponent extends to greater radii (Truiillo fc Browrl 2001). 
Thus the Sun's classification as a narrow-disk system is 
tentative. 

5. SUMMARY 

We present the first optical scattered light images of 
debris disks surrounding relatively old main sequence F 
and K stars. Material around HD 139664 is concentrated 
at 83 AU radius, with a distinct outer edge at 109 AU, 
and a depleted, but not empty, region at <83 AU radius. 
Dust surrounding HD 53143 has a monotonic q~ 1 varia- 
tion in grain number density, and the disk edges from 55 
AU to 110 AU are sensitivity-limited values. The differ- 
ent radial widths appear consistent with a more general 
grouping of debris disks into either narrow or wide ar- 
chitectures. These two categories are probably an over- 
simplification of significant diversity in the formation and 
evolution of debris disks. Future observations should test 
for common traits among these stars, such as stellar mul- 
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tiplicity and the existence of planets. 
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Fig. 1. — False-color, ACS/HRC F606W (A c = 591 nm, AA = 234 nm) images of scattered light from debris disks surrounding (a) HD 
53143 and (b) HD 139664. North is up, east is left and the circular black mask in each image has radius 3". The periphery of each field also 
shows the 3.0" diameter circular occulting mask located at the tip of an occulting bar. Cumulative integration times are 2340 s and 2184 s 
for HD 53143 and HD 139664, respectively. HD 53143 and HD 139664 were observed on 2004 September 11 and October 14, respectively. 
The stellar full-width at half-maximum is 63 mas. The images are corrected for distortion and have 25 mas per pixel. 




Fig. 2. — (a) Mirror averaged image of HD 139664. The white bar has length 1" and the right tip is located at 3" radius from the 
star (the radius where we begin making measurements of disk properties), (b) The same data after a model disk has been subtracted. 
Measurements of the PSF residual along the disk midplane are now indistinguisheable from all other position angles and we consider 
the model a satisfactory fit to the data. Model disks with inclination range i = 85° — 90° fit the observed morphology and midplane 
surface brightness. This range of inclinations is allowed because a vertically fat, edge-on (i = 90°) disk can have isophote morphology that 
resembles a vertically thin disk with i = 85° . (c) Illustration of the narrow-belt architecture by inclining the model HD 139664 disk to 
face-on. Model properties are discussed in Fig. 3. The dust-free hole within 60 AU must be confirmed by future observations. Summing the 
light between 82 and 84 AU radius in our face-on model disk, we obtain m v = 20.3 mag. The perpendicular optical depth is approximately 
10 ((4.6-20.3)/-2.5) = 5 2 x 10 -7 for albedo = 1 
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Projected Radius (arcseconds) 

Fig. 3. — Radial surface brightness distribution along the HD 53143 and HD 139664 disk midplanes. Filled circles and solid lines 
trace the observed disk midplanes, and open circles and dashed lines trace model dust disks inserted into the data orthogonal to the real 
midplane direction. The real and model disk midplanes were sampled in a strip 1.25" wide for HD 53143, and 0.25" wide for HD 139664. 
The noise is indicated by the scatter of data points. A turnover in surface brightness between 4.5" and 5.0" (79 - 88 AU) is common to 
all PSF subtractions of HD 139664 (solid blue circles). If the HD 139664 disk extended radially inward, then the surface brightness profile 
should be continuously brighter inward. In order to match the midplane profile, the model disk (black open circles) has no dust within 
3.40" radius (59.5 AU). Between 3.40" and 4.75" (59.5 - 83.1 AU) the dust density increases as q+ 3 . Between 4.75" and 6.25" (83.1 
- 109.4 AU), the dust number density decreases as q~ 2 ' 5 . The model disk has no dust beyond 109.4 AU. For comparison, a model disk 
(green open squares) with outer radius twice that of the best fit case would be detected as far as 10" radius. The models for HD 54143, on 
the other hand, do not indicate significant variations in radial struct ure between the detect ion limits of 3"-6". All models assume isotropic 
scattering with scaling parameters similar to those used for /3 Pic by Kalas & JcwitJ I1995T , 
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TABLE 1 

Star (rows 1-7) and disk (rows 8-15) properties 





Tin C Q 1 a Q 


TT"Pl 1 QOfiR/l 
XliJ 1030041: 


^ge V^J 1 / 


1.0±0.2 


3 +0 - 7 
u.o_ 2 


Spectral Type 


Al V 


r 01 V - V 


lviass ^ivi0 j 


U.o 


1 Q 
l.o 


T „ IK) 


5224 


6653 


Luminosity (Lq) 


0.7 


3.3 


Distance (pc) 


18.4 


17.5 


m v (mag) 


6.30 


4.64 


Peak disk surf, bright, (mag/arcsec 2 ) 


22.0±0.3 


20.5±0.3 


Disk position angle (degrees) 


147±2 


77 ± 0.5 


Inclination (degrees) 


40 - 50 


85 - 90 


Disk number density gradients (q a ) 


-1 


+3.0 & -2.5 


Inner dust depletion (AU) 


< 55 


83 


Maximum outer radius (AU) 


> 110 


109 


Optical depth from IRAS data a 


2.5 x 10" 4 


0.9 x 10~ 4 


Optical depth from HST data b 


> 1.6 x 10" 5 


1.0 x 10" 5 


Total Dust Mass (g) c 


> 7.1 x 10 23 


5.2 x 10 23 



a The fractional dust luminosity iZuckerman Sonel2004D . 

^Derived from the model disks. Since the disks are optically thin, we sum the cumulative light from the model disk, m^, and quote optical depth 

assume albedo— 1 and the optical depth will scale inversely with the assumed albedo. 
c Dust mass follows from the HST optical depth and assumes a uniform particle radius of 30 /im, density 2.5 g cm -3 and albedo— 1.0. 
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TABLE 2 

Debris disk architectures from scattered light properties' 1 



Name 


Tin (AU) 


Tout (AU) 


Width (AU) 


d (pc) 


SpT 


References 


/3 Pic 


-90 


>1835 


>1745 


19.3 


A5V 


1, 2 


HD 32297 


<40 


>1680 


>1640 


112 


AO 


3, 4 


AU Mic 


~12 


>210 


>198 


9.9 


MlVe 


5, 6 


HD 53143 


<55 


>110 


>55 


18.4 


K2V 


This work. 


HD 107146 


~130 


>185 


>55 


28.4 


G2V 


7 


HD 139664 


83 


109 


26 


18.5 


F5V 


This work. 


Fomalhaut 


133 


158 


25 


7.7 


A3V 


8 


HR 4796A 


60 


80 


20 


67.1 


A0V 


9 


Sun b 


25 


50 


25 




G2V 


10, 11 



Not e. — REFERENCES: (l ) IPantin et alj (TM%: (2) ILarwood &; Kalajl 120011) ; (3) ISchneid er et alj I 2005B: (41 iKalad I2005bfl: C5i lKrist et alj 
I2005D: (6) IKalas et"7fl 120041); (7) lArdila et alj i200ll : (8) IKalas et aTT^OOSall : (9) ISchneider et all 119991) : (10) ILiou & ZooU ll999D : (111 
IMoro-Martin & Malhotral I 120*021) 



a Wc do not include systems younger than 10 Myr that arc likely to possess significant primordial circumstellar gas. Column 2 gives the inner disk 
radius corresponding to the approximate peak of dust number density. Column 3 gives the outer radius. 
From simulations of relatively large grains trapped in resonances with Neptune. 



